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ABSTRACT 

' Non-axisymmetric forces are presented for a sample of 107 spiral galaxies, of which 31 are barred (SB) and 53 show nuclear 
, activity. As a database we use JHK images from the 2 Micron All Sky Survey, and the non-axisymmetries are characterized by 

■ the ratio of the tangential force to the mean axisymmetric radial force field, following Buta & Block (2001). Bar strengths have 
O ^. an important role in many extra-galactic problems and therefore it is important to verify that the different numerical methods 
D ' applied for calculating the forces give mutually consistent results. We apply both direct Cartesian integration and a polar grid 

, integration utilizing a limited number of azimuthal Fourier components of density. We find that bar strength is independent 
^-H ' of the method used to evaluate the gravitational potential. However, because of the distance dependent smoothing by Fourier 
, decomposition the polar method is more suitable for weak and noisy images. The largest source of uncertainty in the derived 
' bar strength appears to be the uncertainty in the vertical scale-height, which is difficult to measure directly for most galaxies. 
' ^ ] On the other hand, the derived bar stength is rather insensitive to the possible gradient in the vertical scale-height of the disk 

■ or to the exact model of the vertical density distribution, provided that a same effective vertical dispersion is assumed in all 
models. In comparison to the pioneering study by Buta & Block (2001), bar strength estimate is here improved by taking into 

' account the dependence of the vertical scale-height on the Hubble type: we find that for thin disks bar strengths are stronger 

than for thick disks by an amount which may correspond to even one bar strength class. 
' We confirm the previous result by Buta & Block (2001) and Block et al. (2001) showing that the dispersion in bar 
, strength is large among all the de Vaucouleurs's optical bar classes. In the near-IR 40 % of the galaxies in our sample have 
bars (showing constant phases in the m=2 Fourier amplitudes in the bar region), while in the optical 1/3 of these bars are 
obscured by dust. Significant non-axisymmetric forces can be induced also by the spiral arms, generally in the outer parts of 
the galactic disks, which may have important implications on galaxy evolution. Possible biases of the selected sample are also 
O ' studied: we find that the number of identified bars rapidly drops when the inclination of the galactic disk is larger than 50". 
-j-J ' A similar bias is found in The Third Reference Catalogue of Bright Galaxies, which might be of interest when comparing bar 
^ \ frequencies at high and low redshifts. 

^ [ Key words: galaxies: spiral ~ galaxies: active - galaxies: statistics - galaxies: kine- 

•'^ ■ matics and dynamics 
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1 INTRODUCTION 

Bars consist mostly of old stellar population (de Vaucouleurs 1955; Elmegreen & Elmegreen 1985), which stresses their 
significance as dynamically important components in galaxies. In fact, a large fraction of galaxies have bars (Block & Wainscoat 
1991; Knapen et al. 2000; Eskridge et al. 2000, Block et al. 2001), indicating that they must be long-lived phenomena in 
galaxies. Bars are fundamental in galaxy evolution, suggested to be driving forces for star formation, forniatiou of riug.s and 
global spiral density waves, and even for the onset of nuclear activity. When quantified the correlations between bar strength 
and the other properties of the galaxies can be studied. The wavelength that best traces the dynamical mass is the near-IR, 
where the obscuration of dust is also less significant than in the optical region. For example, galaxies like NGC 5195, which 
are irregular in the optical may have regular grand-design spiral arms in the near-IR (Block ct al. 1994), which emphasizes 
the importance of a new, more dynamical picture of the morphological structure in galaxies. A step toward that direction is 
the new dust penetrated classification of spiral arms in the near-IR (Block & Puerari 1999; Buta & Block 2001; Block et al. 
2001), in which bar strength plays an important role. 

As discussed by Buta & Block (2001; BB from hereon) there are many quantitative parameters which can be used to 
estimate bar strengths, such as bar-interarm contrast (Elmegreen & Elmegreen 1985) or light remaining after the disk and 
bulge components are subtracted (Seigar & James 1998). The most commonly used method is the maximum ellipticity of a 
bar, an approach justified by the analytical models by Athanassoula (1992), who showed that the non-axisymmetric forces 
in the bar correlate with the bar ellipticity. This method has been recently refined by Abraham & Mcrrifiold (2000), who 
consider both the inner and outer contours of the image to better resolve the ellipticity of a bar. However, the ellipticity is 
not a full description of bar strength. In fact, a more physical approach has been taken by BB who estimate bar torques by 
calculating tangential forces in the bar region, taking into account also the underlying axisymmetric potential. Indeed, when 
refined, the bar torque method is probably the most promising way of estimating bar strengths. 

When the bar torque method is finely tuned, future refinements will include the complex bar structures seen in many 
galaxies; taking more properly into account the vertical scale heights and their gradients, as well as taking into account bulge 
stretch scenarios upon deprojection of the images. In the bar torque method there are also different ways of evaluating the 
gravitational potential and it is important to veriiy that the different methods give mutually consistent results. For example, 
BB used the 2D Cartesian integration method by Quillon ct al. (1994), whose new contribution in the potential evaluation was 
that the vertical density profile of the disk was taken into account in the convolution function. The potential was calculated 
in Cartesian grid by applying Fast Fourier Transform techniques (see also Elmegreen et al. 1989). On the other hand, in our 
study of IC 4214 (Salo ct al. 1999) we evaluated the barred potential by first "smoothing" the image by calculating the Fourier 
decomposition of the surface density in a polar grid. In principle these two methods should give similar results. 

In the present study bar strengths are calculated in JHK-bands for 107 spiral galaxies using the polar method (Salo 
et al. 1999) for the evaluation of the gravitational potential. The method has been improved by taking into account the 
recent observational work showing that bars in early-type galaxies are thicker than in late-type galaxies. Also, the effects of 
a distance dependent scale height, detected in many boxy/peanut shaped disks, are estimated. The algorithm of calculating 
forces is described, the different ways of estimating the gravitational potential are compared and Fourier analysis is applied 
for the analysis of bars. Also, bigises of the sample are studied and the distributions of bar torques among the de Vaucouleurs's 
optical bar classes are compared. In future the method will be further developed to better take into account the observational 
properties of bars and bulges in galaxies. The measurements of this paper have been used in comparisons of bar strengths 
between active and non-active galaxies by Laurikainen et al. (2002). 

2 THE SAMPLE SELECTION 

A sample of spiral galaxies was selected from the Third Reference Catalogue of Bright Galaxies (de Vaucouleurs et al. 1991, 
RC3) requiring that Bt < 12.5 mag, cz < 2500 km and the inclination INC < 67°. As a database we use the 2 Micron All 
Sky Survey (2MASS). As only about 50 % of the galaxies in our original sample have images available in 2MASS, the sample 
in the present form is not magnitude-limited. Some of the weakest images were also eliminated, because bar strengths could 
not be measured in a reliable manner for them. The final sample consists of 107 spiral galaxies, of which 31 are barred (SB) 
in RC3, 53 show nuclear activity, 42 are early-type spirals (SO/a-Sb) and the rest of the galaxies belong to late- types. To 
active galaxies we include Seyferts, LINERs and Hll/Starburst galaxies, the type of activity being taken from NASA/IPAC 
Extra-galactic Database (NED), where the most recent classifications are available. 

The frequency of bars in our sample is typical in comparison to other samples of galaxies. Including to barred galaxies 
both SB and SAB types we find a bar frequency of 62 %, which is similar to that obtained by Sellwood & Wilkinson (1993) 
for field galaxies (60 % barred) and by Ho et al. (1997) for a magnitude-limited sample {Bt <12.5) of spiral galaxies (59 % 
barred). A somewhat larger fraction of barred galaxies has been found by Hunt & Malkan (1999) for a sample selected by 
12 nm radiation (69 % barred) and by Moles et al. (1995) for a magnitude-limited sample extending to Bt=13 mag (68 % 
barred). The numbers of SB and SAB galaxies in our sample are rather similar (29 % and 33 %). While identifying bars in 
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the near-IR by Fourier techniques (sec Section 5.3) the fraction of barred galaxies was found to be 40 %. 

The fraction of SB galaxies among Seyferts and LINERs in our sample is 30 %, being similar as for the whole sample, in 
agreement with the bar fractions detected in the samples by Ho et al. (1997) and Mulchaey & Regan (1997). An enhanced 
frequency of bars is generally associated to Hll/Starburst galaxies, for example, the fraction of SB-galaxies among the 
Markaxian Starburst galaxies is even 75 % (Hunt & Malkan 1999). The bar fraction among HH/Starburst galajcies somewhat 
decreases when weaker nuclear star formation activity is considered: namely, in the whole 12 /nm sample by Hunt & Malkan 
the fraction of SB-galaxies is 53 %, but while concentrating to smaller distances (cz < 5000 km s^^), it decreases to 46 %. In 
our sample only 2 out of 17 HH-galaxies are classified as SB, which can be partly understood by the small distances of the 
galaxies (< cz >= 1200 km s~^): most probably we are picking the lower end of the nuclear HH-luminosity function, where 
bars may not be the driving forces for nuclear star formation. 

The frequency of active nuclei in galaxy samples depends strongly on the redshift range studied: at small distances almost 
all gala^xies show nuclear activity in some level (Ho et al. 1997), while the number of strong active nuclei increases with redshift. 
In our sample 53 % of the galaxies have active nuclei in terms of Seyferts, LINERs and Hll/Starburst galaxies. As we have the 
same apparent magnitude limit as Ho et al., more Hll/Starburst galaxies would be expected. The reason why we don't have 
more Hll/Starburst galaxies most probably is that NED may not sample well the lowest level activity. The morphological 
types of Seyferts and LINERs in our sample are peaked to Sab galaxies, which is between the mean morphological types for 
Seyfert 1 (Sa) and Seyfert 2 (Sb) gala:xies (Malkan et al. 1998), whereEis 3/4 of the non-active galaxies belong to late Hubble 
types. 



3 EVALUATION OF THE GRAVITATIONAL POTENTIAL 

Like in BB, also here the bar strengths are estimated from the magnitude of the non-axisymmctric gravitational perturbation 
in comparison to the mean axisymmetric radial force field. For each radius r and azimuth we calculate the tangential force 
Ft = ^d^/d4> and the radial force Fr = d^/dr, and define the relative strength of the perturbation as 

Qr{r) = FT""'''{r)/ < FR{r) >, (1) 

where the average of radial force over azimuth is taken. Following BB the tangential force maxima were calculated separately 

in four quadrants of the image and the mean of these in each distance is used as the maximum Ft- In order to obtain a single 
measure for the strength we use Qb, which is the maximum Qt in the bar region. The radial distance where this mgiximal 
perturbation takes place is denoted by rgb- 

For the force calculation the gravitational potential $ in the bar region must be evaluated. We assume that the surface 
density E is proportional to the surface brightness obtained from the near-IR image, and that the vertical density distribution 
follows a model profile Qz{z), normalized to unity when integrated over z. First it was assumed that the z-dependence of the 
mass density obeys the same formula everywhere in the galaxy (see Section 4.1), in which case the gravitational potential in 
the central plane of the galaxy can be written as (e.g. Quillen et al. 1994) 



/ ^{x ,y')g{x - X ,y - y')dx'dy' 

-OO J — OO 



(2) 



where the integral over z-direction is included to the convolution function g{x — x',y — y') = g'(Ar'), with Ar^ = (a; — x'Y + 
{y — y'Y , defined as 



/OO 
Q.(z){/^r'' + z^)-^'^dz. 
-OO 



(3) 



Several models for (z) were applied, including the exponential model with scale-height hz 



Q-Az) = i^J'xp{-\z/hz\), (4) 
and the often used isothermal sheet model (van der Kruit & Searle 1981), 

= 7;r- — sech'^{z/hsech2}, (5) 

where hsech2 is the isothermal scale-height. Usually, when isothermal model is used, hsech2 is set equal to 2hz, to yield a same 
slope at large z as for the exponential model. According to van der Kruit (1988) galactic disks may, however, deviate from 
isothermal shape near the galactic plane, in which case a better approximation for the density distribution would be 

Qz{z) = —r^ Sech{z / hsech) ■ (6) 

^ '^sech 

In the case of sec/i-law, Barnaby & Thronson (1992) identify hsech with {n/4)hz, to get the same central plane density as in 
the sech^ -law. However, in this study we briefly check sech-law with hsech = hz, in which case the slope of the vertical profile 
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corresponds to exponential model at large z. De Grijs et al. (1997) have explored a more general family of fitting functions, 
which include the above functions as special cases. As will be discussed below the adopted vertical model has rather little 
influence on the derived forces, provided that appropriately defined scale-heights are used. 

In comparison to BB we use a somewhat different method for the evaluation of the gravitational potential, mainly because 
the 2MASS images used here have a more poor resolution (l"/pix) and lower signal-to-noise (S/N) ratio than the images used 
by BB. Thus, instead of calculating the gravitational potential directly from the image pixels, the images are first "smoothed" 
by calculating the Fourier decompositions of the surface densities in different radial zones, 



m— oo 



<f>) = Eo(r) ^ E„(r, = Ao(r) + ^ Am{r) cos [m (0 - 4>m,(r))] . (7) 

m=l m=l 

The density of each Fourier component is then separately tabulated as a function of radius and azimuth, and for each density 
component the corresponding potential component is obtained by 

/•oo r-2TT 

^m{r,<t>,z = 0) = -G r'dr' S^(r', 0')5(Ar)d0', (8) 



with Ar^ = r''^ +r'^ — 2rr'cos{(j)' — (j))- The integration over azimuthal direction is carried out by Fast Fourier Transform (FFT), 
whereas in the radial direction a direct summation is used. An azimuthal offset of one half bin is used between density and 
potential locations, and the force components at density locations are obtained from the potential by numerical differentiation. 
For test purposes, we also applied Cartesian potential evaluation, solving Eq. (1) with 2D FFT as in BB. In comparison to 
direct Cartesian force evaluation the advantage of our polar method is that also rather weak and noisy structures can be 
measured in a reliable manner, due to the smoothing implied by Fourier decomposition. Also, it is possible to limit the density 
to even components, most likely to characterize the non-axisymmctry related to the bar. Our method also gives directly the 
different Fourier modes of the potential and force components, which are sometimes of interest. 

Like in BB and Quillen et al. (1994), we made use of the fact that the convolution function g{Ar) can be numerically 
integrated and tabulated as a function of Ar/h, h denoting the vertical scale- factor used, so that the integration over z- 
direction can be replaced by a much faster interpolation from prc-calculated tables. In addition, since in the polar method 
only the integration over azimuth is carried out by FFT, a distance-dependent h can be used: in this case g{Ar/h) in Eq. (8) 
remains cyclic with (j>' — 4> even if /i = h{r'). In the Cartesian 2D FFT this can not be done as the convolution function needs 
to be cyclic in both x — x' and y — y' . 

In Fig. 1 we display the convolution functions corresponding to various vertical models. The upper row shows the density 
profiles and convolution functions corresponding to exponential, sech, and sec/i^-model profiles. The difference in g{Ar) is 
significant only for Ar < h^.ln the lower row of the figure we illustrate the relative insensitiveness of g{Ar) for a wider range 
of Qz{z) models, including the Gaussian model, 

^'^^"^ = ~7t^ ^ ) ' 

y ^T^n gauss \ ^'f^gauss J 

and a uniform slab-model, 

ez{z) = , -hunt < Z < huni- (10) 

In order to make the models comparable, the scale factors in each case were chosen in a manner that yields the same vertical 
dispersion as the exponential model with h = hz, namely < z^gz{z) > / < Qz{z) >= 2hz^ . In this case 

hgauss I hz — - 

hsech2/hz = V24/7r, 

hsech/hz = VS/tT, (11) 
huni/hz = a/6. 

As expected, the more peaked the density profile is toward the central plane, the larger is the value of g{Ar) when 2: — > 0. 
In each case g{Ar) oc 1/Ar for large Ar, while g{Ar) increases logarithmically when Ar — > 0. In the case of a uniform slab, 
analytical integration yields 



log ('""' + ^5-^1, (12) 
huni \ Ar / 
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in agreement with the result of the numerical integration shown in Fig. 1. For comparison, Fig. 1 displays also the case of 
using 2-D softcucd gravity, 



In conclusion, the convolution function depends very little on the model used for the vertical mass distribution as long 
as models with a same vertical dispersion are compared. In the next section bar strengths calculated using different vertical 
models are compared, as well as results obtained by the two different integration methods. 

4 TESTING THE ALGORITHM 
4.1 Polar vs. Cartesian integration 

Wo next investigate our algorithms of force evaluation using NGC 1433 as a test case. For this galaxy Ron Buta has provided 
us with his high quality H-band image, with a 1.141 arcsec pixel scale (Buta et al. 2001). Fig. 2 compares our standard polar 
method, and our Cartesian potential evaluation, where the de-projected image is interpolated to a density array and then FFT 
in Cartesian coordinates is applied with a grid resolution ranging from 64 x 64 to 512 x 512. In both polar and Cartesian cases 
the forces where calculated from a region with radius of rmax = 180 pixels (205 arcsecs). According to Fig. 2 the obtained 
Qr-profiles converge rapidly as more Fourier components are included, or when a finer Cartesian grid is used. The profiles 
obtained with the two methods also agree well: the maximum Qt agrees within a few percent. The only differences appear at 
large r, where the polar method yields smaller Qt- This is due to distance-dependent smoothing implied by the polar method, 
whereas in the Cartesian integration the maximum tangential forces for large r are mainly due to spurious values connected 
to individual image pixels. Thus, since the outer parts of the galaxy have a low S/N-ratio, the smoothing in the polar method 
is also physically well motivated. 

As the polar method is based on calculating Fourier amplitudes for different density components, it is interesting to check 
how much each component contributes to the total force: in Fig. 2 even components up to m=10 are studied. We have used 
an exponential vertical profile, both with hr/hz = 2.5, representing an early-type gala^xy, and with hr/hz = 10, to study the 
effect of assuming a very flat disk. The density amplitudes m=2 and m=4 are known to be the strongest in bars, while in 
early-type galaxies the m=6 and m=8 may also be significant (Ohta 1996). For NGC 1433 (type Sb) the density amplitudes 
we find in H-band are fairly similar to those displayed for the I-band by Buta (1986), with A2/A0 attaining a maximum 
0.97 at r = 110", and a secondary maximum of w 0.6 at r w 50". For m = 4,6,8,10 the maximum Am/Ao are about 
0.42, 0.22, 0.13, 0.09, respectively, and for the resulting force the maximum appears at rqb = 70", which corresponds to about 
90% of the bar radius {Rbar = 77", Buta 1986). For hr/hz=2.5 about 75% of the tangential force is due to the m=2 density 
component, and with the inclusion of the m=4 component Qt increases to 97 % of its value obtained by including all even 
components up to m=10. However, for the case of a flatter galaxy the influence of higher order Fourier components would 
become somewhat more important: in the case of hr/hz = 10, Qb = 70%, 90%, 98% of its maximum value, when including 
the density components up to rrimax = 2, 4, 6, respectively. Thus even in this case the density components above m=6 have 
fairly little influence. This difference in the relative importance of various Fourier components follows from the fact that 
planar density variations correspond to force variations only if their planar scale significantly exceeds hz'- smaller hr/hz thus 
suppresses the force variations corresponding to large m. 

In the above example the calculation region covered well the optical disk of the galaxy, and the maximum Qt was 
obtained at rmax/rqb ~ 3. However, in the case of 2MASS images the outer disks for distant galaxies are sometimes not deep 
enough. In order to check the influence of the size of calculation region, a series of integrations with decreasing rmax was 
conducted for NCC 1433. For rmax/rqb = 2 the resulting value of Qb was only about 2% smaller (no shift in rgt), whereas 
for rmax/rqb = 1-5 the reduction was already significant, amounting to 15%. About 2/3 of this reduction is due to reduced 
tangential force, while the remaining part comes from the overestimated radial force due to disk truncation. However, the 
location of Qr-maximum was still fairly little affected (66" vs. 70"). All our subsequent force evaluations from 2MASS images 
are based on images covering a radius at least twice as large as the derived rqt- 

We also checked the influence of different model functions for the vertical profile on NGC 1433 (Fig. 3). The models 
correspond to those studied in Fig. 1, each having the same vertical dispersion. As expected, based on the behavior of the 
convolution functions, a more centrally peaked vertical profile yields a slightly larger Qb, although the shape of the Qr-profile 
is little affected. For early-type galaxies with hr/hz = 2.5, the difference in Qb between the exponential and uniform models is 
15%, which is much less than the difference between the cases hr/hz = 2.5 and hr/hz = 10 for a fixed vertical model function. 
Also, while comparing more realistic vertical models, namely the exponential and the isothermal model, Qb is affected only 
5%. In general, when the vertical extent of the disk is reduced, the exact form of the vertical density law becomes less 
important. In the limit of very large hr/hz the difference between including 3D vertical profile and the use of softening in 2D 
force-evaluation becomes also small, in agreement with Salo et al. (1999). However, for large disk thicknesses {hr/hz = 2.5), 
the use of softened gravity severely underestimates Qb- 
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4.2 Distant dependent disk thickness 

As many galaxies have boxy/peanut shaped structures, for whieh the vertical thickness has found to increase towards the 
outer parts of the galaxies, we also studied a case with distant dependent vertical scale-height. Radial distributions of Qt 
were calculated for NGC 1433 using an exponential density model, when hz was assumed to have a mean value and a gradient 
typical for early-type galaxies (see Section 5.1). As shown in Fig. 4 the vertical density gradient has a fairly small effect on 
bar strength: in comparison to the case where a constant hz was assumed, a gradient of 0.05 in hz affects Qt at maximum 
only by 4%. In the figure three cases are demonstrated showing how the location of the mean vertical thickness aflFects Qt 
measurement. In comparison to constant vertical scale-height, radially increasing hz typically slightly reduces the derived bar 
strength. This is due to two competitive effects: on one hand the tangential force increases in the region where the vertical 
scale- height is reduced, but at the same time axisymmetric radial force is also increased in the same region (see Fig. 5). 

In interacting galaxies gradients of the vertical scale-heights are generally larger than for normal galaxies, and contrary 
to normal galaxies, may also appear in late-type systems. Therefore we studied also some cases, where the hz gradient was 
twice as large as for normal galaxies. For hr/hz=2.5 the eflfect of the gradient on Qh was still unimportant Qh being reduced 
only by 3%. Thin disks were even more less affected: for hr/hz = 10 we found that while adding a gradient of 0.03, Qt was 
maintained unchanged. We can safely conclude that the vertical scale-height gradient, typical for boxy/peanut shaped bars, 
does not significantly affect bar strength. 

4.3 Comparison to BB 

We next compare our results with those obtained by BB. Bar strengths for 13 galaxies in the sample by BB were calculated in 
H-band, using the same orientation parameters as used in BB (see Table 4). Also, following BB the vertical mass distribution 
was approximated by an exponential function with hz —325 kpc, which is the scale-height of the disk in the Milky Way. As 
shown in Fig. 6, there is a small shift toward somewhat larger Qt in our measurements, but generally the agreement between 
the two measurements is quite good. No softening was used in our potential evaluation. A possible cause of the small difference 
is that some form of additional gravity softening was included to the potential evaluation in BB, besides using a vertically 
extended density model. This possibility is suggested by the form of the convolution function g{r) shown in Fig. 2 in Quillen 
et al. (1994), being the method used by BB. In comparison to our g{r) (Fig. 1), their function seems to attain a constant 
value for Ar/h -> 0, much like our curve for softened gravity. Nevertheless, no explicit gravity softening was mentioned in 
Quillen et al., so it is not clear whether this is the cause for the difference. It is also worth noticing that the comparison was 
made using 2MASS images, many of which were of more poor quality than the images used by BB. In Table 4 we also show 
that the vertical density model barely affects the derived bar strength. 

5 CALCULATION OF THE NON-AXISYMMETRIC FORCES FOR THE SAMPLE 
5.1 Observed parameters in the potential evaluation 

As discussed in the previous section, several approximations were made while calculating the gravitational potential, largely 
following BB. The main assumption was that the mass-to-luminosity ratio (M/L) is constant throughout the disk. This 
assumption was made for simplicity, because otherwise possible M/L- variations along the disk should be known for each 
individual galaxy. The studies of color gradients in galaxies have shown that the central regions of the disks are often redder 
than the outer disks, which actually suggests larger stellar M/L-ratios in the inner disks (Bell & de Jong 2001). Another 
indicator of a possible non-constancy of the M/L-ratio comes from the comparison of the optical surface photometry with the 
surface mass densities obtained from the rotation curves, showing that the M/L-ratio may vary along the disk (Takamiya & 
Sofue 2000). However, it is possible that the M/L-ratio is rather constant in the bar region (see Quillen et al. 1994), in which 
case the made assumption should be reasonably good. 

In the evaluation of the barred potential a model for the vertical mass distribution was assumed. The most commonly 
used models are the exponential and isothermal functions, both being physically justified. An isothermal density function is 
expected if stars, once formed, do not interact with the other components of the galaxy. Stellar populations of different ages 
can then be understood as quasi-independent components with different velocity dispersions (Dove & Thronson, 1993) so that 
a more complete picture would be achieved by assuming a super-position of a large number of isothermal sheets with different 
hz (Kujken 1991). On the other hand, if gas settles into an equilibrium prior to star formation or considerable heating of the 
disk has occurred during the life of a galaxy, an exponential density profile is expected (Burkhead & Joshii 1996). However, 
empirically the two functions are difficult to distinguish, because the density functions look similar at large vertical heights, 
whereas near to the galactic plane the evaluation of the density function is complicated by the effects of dust. For the vertical 
density distribution we used an exponential function. However, based on discussion in Section 4.1, the uncertainty in the 
vertical model function should have fairly small influence on derived forces. 

For bar strength measurements, more critical is the thickness of the disk. It has been found to depend on the morphological 
type of the galaxy (de Grijs 1998) and also to correlate with the radial scale-length of the disk (Wainscoat et al. 1989; van 
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der Kruit 1988). For early-type galaxies we used hr /hz=2.5, and for late-type systems /ir//iz=4.5. For hr wc use the optical 
V-band scale-lengths from Baggett et al. (1998), and if not available, they were estimated from the 2MASS images by us. 
Optical scale-lengths were used, but they are expected to be rather similar with those in the near-IR: namely using the bulge- 
to-disk decompositions given by de Jong (1996) for 186 spiral galaxies, we found that the scale-lengths in V-band deviate on 
the average only 5% from those in K-band. For three of the galaxies, pgc 10266, pgc 15821 and pgc 40097, the scale-lengths 
by Baggett et al. were judged unrealistic: for two of the galaxies they were rather measures of the brightness slopes in the 
bulge region, and for one galaxy the given scale-length represented the outermost very shallow part of the disk, while wc are 
interested in the disk under the bar. Therefore, also for these three galaxies the scale-lengths were estimated from 2MASS 
images. 

Bars and bulges have often rather complex structures and may in some cases be difficult to distinguish from each other. 
For example, 45 % of all galaxies (SO-Sd) may have boxy or peanut shaped bulges or bars, having vertical disk profiles that 
become thicker towards the outer parts of the disks (Lutticke et al. 2000, Schwarzkopf & Dettmar 2001). These gradients 
arc pronounced in early-type galaxies, but generally do not appear in latc-typc galaxies. The boxy/peanut structures are 
often thought to be bulges, but resent observations rather support the idea that they more likely represent the thick parts of 
the bars. This interpretation is also supported by the simulations by Athanassoula (2002) who has shown that boxy/peanut 
structures are formed from the particles of the disks or bars during the evolution of the galaxy. In this work we have made 
the simplified assumption that hz is constant throughout the disk, although the effect of the thickening of the disk toward 
the outer parts of the disk is also investigated. In the test we used a dhz/dr gradient of 0.05, typical for normal SO/a galaxies 
(Schwarzkopf & Dettmar 2001). Galaxy interactions and mergers of small satellite galaxies can also be efficient in thickening 
the disk especially in the outer parts of the galaxies (Toth & Ostriker 1992; Walker et al. 1996). In fact, observations by 
Schwarzkopf & Dettmar have shown that the vertical heights of the disks for interacting galaxies can be twice as large as for 
non-interacting galaxies. Bulges were assumed to be as flat as the disks, which is also a simplified assumption, but may still 
be valid for the triaxial bulges of SB-galaxies (Kormendy 1982, 1993). Evidently, the treatment of bulges need to be improved 
in future. 

5.2 Calculation of forces for sample galaxies 

As a database JHK images of the 2MASS survey were used. The spatial resolution of the images was one arc second, and the 
image quality was generally best in the H-band. Our procedure of estimating the non-axisymmetric forces consisted of the 
following steps: (1) cleaned sky-subtracted mosaics were constructed, (2) galaxies were de-projected to face-on orientation, 
(3) the images wore robinnod by a factor of two, (4) Fourier decomposition of the surface density was calculated, and barred 
potentials were evaluated using the even components up to m=6, (5) the tangential {Ft) and the mean ajcisymmetric radial 
forces (< Fr >) were calculated, and finally (6) maps of the force ratios were constructed: 

QrirA) =FT{r,(l,) / < FR{r,,p) >. 
For a bar the map shows four well defined regions where the force ratio reaches a maximum or minimum around or near the 
end of the bar. As in BB, we call this structure as a "butterfly pattern" . 

The field of view in the 2MASS images is relatively small and also a large fraction of the galaxies in our sample arc quite 
nearby objects so that mosaics of 2-5 images were generally made. The image quality was not as good in the borders as in 
the central parts of the frames, which in principle could seriously affect the quality of the mosaics in the regions of interest. 
For the most nearby galaxies this was not a problem, because the S/N-ratio was high in all parts of the images. For the 
distant objects the problem was solved so that we never combined images if their borders appeared in the bar region or in 
the central parts of the galaxies. Before combining the images the over-scan regions were removed and the sky values given 
in the image headers were subtracted. The background levels of the frames were then refined to get similar count levels in 
the galaxy regions in the combined images. Positioning of the frames in the mosaics was done using stars common in the 
combined images. Finally the foreground stars and bad pixels were rejected. 

The cleaned (mostly mosaics) images were de-projected to face-on orientation using the position angles (PA) and incli- 
nations of the disks given in Table 1. In the table some other properties of the galaxies like the mean revised morphological 
type Tm, the blue apparent magnitude Bt, the type of nuclear activity, and the scale- length of the disk hr, are also shown. 
The orientation parameters, the morphological types and the apparent magnitudes are from RC3 if not otherwise mentioned. 
For some of the galaxies, instead of using the orientation parameters from RC3, they were estimated from the Digitized Sky 
Survey Plates by us. For example, for pgc 18258 the surface brightness contours are clearly affected by the super-position of 
a small companion, which has not been taken into account in the orientation parameters in RC3. For the galaxies that had 
no estimation of PA in the literature, but the disks were in nearly face-on orientation, we used INC=0. In order to estimate 
properly the axisymmetric radial forces the bulges were not subtracted. 

The maximum of non-axisymmetric forces, Qt, and their radial distances, rqb, are shown in Tables 2 and 3, measured in 
J, H and K-bands. The uncertainties attached in the tables are the maximum deviations between the four image quadrants. 
However, the largest uncertainty (see Section 4.1) is due to the observed scatter in hr/hz within each morphological type. 
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which for example for Sc-galaxies induces an uncertainty of about 15 % in Qh. Non-ajcisymmetric forces are also sensitive 
to the orientation parameters of the galaxies: BB estimated that an uncertainty of ± 10 % in the inclination and position 
angle can induce an uncertainty of two bar strength classes (one bar strength class corresponds to 0.1 units in Qb). Indeed, 
in future this work can be improved when accurate position angles and inclinations will be available for all of the galaxies 
studied. Especially for interacting galaxies the photometrical orientation parameters are generally determined from waxped 
or distorted outer disks so that the kinematic observations for them give more reliable values. For example, for M51 the 
kinematically (TuUy 1974) and photometrically (Spillar et a. 1992) estimated inclinations deviate by about 20*^, and for IC 
4214 by 10° (Buta et al. 1999), for which the difference can be explained by the bar potential (see Salo et al. 1999). On the 
other hand, the uncertainty caused by the vertical model (exp v.s. sech^) is negligible for thin disks and only about 5 % for 
thick disks. Also, the number of Fourier modes or the size of the measurement region do not affect Qb significantly, if large 
enough measurement regions {vmax > 'irqb) and enough Fourier modes (even modes up to m=6) are used. The effect of the 
boxy /peanut shaped structures in terms of increasing vertical scale- height toward the outer parts of the disks also appeared 
to be insignificant for bar strength. 

5.3 Identification of bars and calculation of bar ellipticity 

Bars in the near-IR were identified by Fourier techniques. In distinction to bars in general, we call "classical bars" the 
morphological structures with the ratio of the Fourier amplitude A2/AQ larger than 0.3, and with the m=2 phase maintained 
nearly constant in the bar region (Table 2). For these bars the m=4 amplitudes arc also pronounced, there is a clear maximum 
in the Qr-profile in the bar region and the "butterfiy pattern" shows four regular structures. The Qr-profiles and the "butterfly 
patterns" for these galaxies are shown in Fig. 7, showing also the m=0,2, and 4 surface brightness profiles. The length of 
the region with a constant m=2 phase was taken to be a measure of the bar length. Both the m=2 and m=4 phases are 
maintained nearly constant in the bar region, but in some cases the m=4 amplitude drops at a slightly shorter distance, in 
which case we used the mean length of the constant m=2 and m=4 phases as the bar length. 

In addition to the "classical bars" , other bar-like structures and non-axisymmetries are identified in many of the galaxies 
in our sample. Actually a large majority of galaxies have non-axisymmetric forces, manifested as maxima in the Qr-profiles, 
whereas in 20 % no non-axisymmetric forces were detected above the background level ((56=0.01-0.1). Due to differences in 
image quality no single minimum Q^-valuc defines the bar-like potential. The non-axisymmetric forces presented in Table 3 
deviate from the "classical bars" in a sense that the m=2 phase is not maintained constant in the assumed bar region and 
the regular "butterfly patterns" do not necessarily appear in the force fleld. In 15 % of the galaxies in our sample the Qt 
maxima are manifestations of strong spiral arms in the outer parts of the galactic disks. 

As discussed by BB, both spherical and flattened bulges can affect Qb-measurements: in a case of an intrinsically spherical 
bulge, the effect of assuming a bulge as thin as the disk is to overestimate the axisymmetric radial force and consequently 
to underestimate the relative bar strength. On the other hand, while de-projecting the image to face-on orientation bulges 
might cause artifacts in the direction of the minor-axis of the disk. The problem of large bulges was here avoided by limiting 
to those cases where the maximum tangential force appeared outside the bulge region. The small bulges still makes it difficult 
to detect mini-bars, but by subtracting the m=0 component and by taking into accomit the dc-projcction effects, mini-bars 
could be detected for pgc 10488, 33371, 37999, 40153 and 43495, being previously identified also by Buta & Crocker (1993), 
Perez-Ramirez et al. (2000), Knapen et al. (1995), and by Block et al. (1994). 

Based on the analytical work by Athanassoula (1992) the maximum ellipticity of a bar can be used as an approximation 
of bar strength, and its radial distance as an estimate of bar length. We calculated the ellipticity profiles using a method 
described in Laurikainen & Salo (2000) in which ellipses were iteratively fitted to the isophotes of the surface brightnesses. 
The maximum ellipticities e and their radial distances in the bar region for the "classical bars" are shown in Table 2 . In 
Laurikainen et al. (2002) these ellipticity measurements were utilized to show a good correlation between Qb and e. 

6 THE SAMPLE BIASES 

The measurements reported in this work are used to compare bar strengths of the active and non-active galaxies (Laurikainen 

et al. 2002), so that it is important to study possible biases between different subgroups in the sample. 

Active galaxies in our sample appeared to be somewhat brighter than the non-active systems, which is illustrated in Fig. 
8a, showing also a weak correlation between the bar length and the absolute blue magnitude Mb, of the galaxy. However, 
while scaling the bar length to the scale- length of the disk, the bias was largely diluted (Fig. 8b). A similar correlation has 
been found previously by Kormendy (1979) for optically measured bar lengths. Due to the magnitude bias longer bars may 
have been selected for the active galaxies in our sample. However, we confirmed that this bias does not affect the mean bar 
strengths in the compared subsamples. This was checked by dividing the non-active galaxies to two magnitude bins with Mb 
larger and smaller than -19.8 mag, resulting to practically identical mean forces, < Qb > = 0.22 ± 0.10 and 0.24 ± 0.12, 
respectively (the uncertainties indicate the sample standard deviations). 
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The relative number of barred galaxies rapidly drops in our sample when the inclination of the disk is larger than 
50° (see Fig. 9). This is the case both for SAB and SB-galaxies and for the "classical bars" identified in the ncar-IR. The 
inclination distribution for the galaxies in the whole sample largely follows a random distribution of orientations, dN/di ~ sin 
i, expected for an unbiased sample. Again, since in Laurikainen et al. we compare bar strengths between different subsamples 
it is important that this inclination bias does not affect the mean <56-values. It appeared that Qt, does not correlate with the 
inclination of the disk (see Fig. 10) : especially for the galaxies with no "classical bars" Qh is similar for all inclinations. For the 
galaxies with "classical bars" there may be a lack of strong bars among the nearly face-on galaxies (INC < 25°) and clearly 
among the highly inclined galaxies (INC > 60°), which is probably a manifestation of the uncertainties in the inclination 
determinations. 

In order to check whether the dependence of bar detection frequency on inclination is specific for our sample or docs appear 
also in larger samples of galaxies, we picked up all galaxies brighter than 15.5 mag from RC3 and constructed histograms for 
SA, SAB and SB-galaxies (Fig. 11), similar to those in Fig. 9. Evidently, a similar bias appears also in RC3: for SA-galaxies 
the inclinations are well sampled to 60°, for SB-galaxies to 50° and for SAB-galaxies the limit of well sampled galaxies is even 
lower than that. Quite surprisingly, the number of SAB-galax;ies at inclinations larger than 40° drops much more rapidly than 
the number of SB-galaxies indicating that classification of a galaxy as SAB is very ambiguous. It also seems that in RC3 there 
is a deficiency of galaxies with very low assigned inclinations (loss than 20"), supporting the above interpretation that the lack 
of strong bars in nearly face-on galaxies is solely due to uncertainties in the orientation parameters. The found inclination 
bias might be of importance for example when bar frequencies are compared between low and high redshift galaxies. In order 
to study a possible redshift dependence of the bias, barred galaxy fractions in two inclination bins were compared at two 
magnitude intervals in RC3. The bias was found to become more significant toward the fainter galaxies: the relative number of 
SB galaxies with INC > 50° dropped from 32 % to 25 % while going from the magnitude interval Bt < 12 mag to Bt= 13.5 
- 14.5 mag, whereas the relative SB galaxy numbers with INC < 50° were identical (38 %) in the two magnitude intervals. 
This means that when bar frequencies at high and low reshifts axe compared, bar frequencies of distant galaxies are easily 
underestimated. 

7 COMPARISON OF OPTICAL AND NEAR-IR BARS 

The distributions of the non-axisymmetric forces among the de Vaucouleurs's (1963) optical classes SA, SAB and SB are 
shown in Fig. 12 in ,1, H and K-bands. In agreement with BB and Block ct al. (2001), the overlap between the different de 
Vaucouleurs's types is significant. Based on estimating the ellipticities of bars, a somewhat different result weis obtained by 
Abraham & Merrifield (2000), who argued that SB-galaxies are clearly separated from SA and SAB-galaxies in bar strength. 
In our sample SA and SAB galaxies have fairly similar non-axisymmetric forces, the only difference being that some SA 
galaxies have no non-cixisymmetric forces. 

Actually, it seems that SB-galaxies can belong to any of the bar strength classes from 1 to 6, as defined by BB. In our 
sample the minimum Qb=0.09 for SB galaxies is largely due to the image quality. We find that bars in SB galaxies are very 
similar to the "classical bars" in the near-IR: namely 95 % of all SB galaxies in our sample are classified as "classical" in 
the near-IR. They also have very similar bar strength distributions in Fig. 12. However, there are 30 % more "classical bars" 
than SB bars, the excess being distributed to all bar strengths, which means that even strong bars can be hidden by dust in 
the optical. In addition to bars significant tangential forces can be induced also by the spiral arms, especially in the outer 
parts of the disks, of which M51 (pgc 47404) is a good example. The spiral related forces may in some cases amount even to 
(56=0.26, corresponding to a bar strength class 3. This is quite interesting, because it means that even spiral arms can induce 
non-axisymmetric forces in the level typically associated with moderately strong bars. These non-axisymmetries may have 
important implications in secular evolution in galaxies, as for example for the onset of near nuclear star formation. 

8 CONCLUSIONS 

Non-axisymmetric forces are calculated for 107 spiral galaxies in J,H and K bands using a method where gravitational 

potentials are evaluated in a polar grid. Non-softened convolution function is applied and the vertical distribution of matter 
is approximated by an exponential function. The vertical scale-height of the disk is taken to be a certain fraction of the radial 
scale- length of the disk, and this ratio is assumed to be larger for the early than for the late-type galaxies. The M/L-ratio is 
assumed to be constant throughout the disk. The vertical mass distribution is generally assumed to obey the same formula 
everywhere in the galaxy, but tests were also performed to estimate the effect of radially non-constant vertical thickness. The 
phases of the Fourier density amplitudes are used to estimate the lengths of the bars. One of the main concerns of this study 
is to verify that the different methods of calculating the gravitational potential give mutually consistent results, most of the 
tests being carried out using a high quality H-image of NGC 1433 (Buta et al. 2001). In comparison to BB our method is 
more suitable for weak and noisy images. Also, it is possible to limit only to even Fourier decompositions, most likely to 
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characterize the non-axi.symmetry related to the bar. Likewise in polar method it is easy to study distant-dependent hz. The 
isophotal ellipticities of bars are also estimated, to facilitate comparisons to bar strengths estimated from maximal forces. 
The main results are the following: 

(1) Cartesian and polar grid methods for the potential evaluation are compared. In the Cartesian method the image 
is sampled to a density array and then 2D FFT in Cartesian coordinates is applied. In the polar grid method Fourier 
decomposition of density is calculated in a polar grid using FFT in azimuth and direct summation over radius. We found that 
similar results are obtained by these two methods for good quality images, provided that enough Fourier components (up to 
m=6) are included, and the resolution of the Cartesian grid is sufficiently large. 

(2) Bar strength is found to be rather insensitive to the vertical mass model of the disk, as long as a same vertical 
dispersion is assumed for all models (e.g. haech2/hz = V^i/n, hsech/hz = \/8/7r). Boxy/peanut shaped structures, in terms of 
non-constant vertical scale-heights along the disk, were also found to be quite unimportant for the evaluation of bar strengths. 
These parameters affect Qb less than 5%. The largest uncertainties in Qb are associated to the large scatter in the observed 
vertical scale-height of the disk within one Hubble type, and to the observed uncertainties in the orientation parameters of 
the disks, which both may induce uncertainties of about 10-15 % in Qb- 

(3) Significant non-axisymmetric forces {Qb >0.05) are detected in 80 % of the galaxies in our sample. In most cases 
they were interpreted as bar-like features, based on significant m=2 Fourier amplitudes in the bar regions and distinct 
"butterfiy patterns" in the Ft/ < Fr >-ratio maps. In 40% of the galaxies "classical bars" were detected, determined as 
having A2/A0 > 0.3 and the m=2 phases maintained nearly constant in the bar region. In some of the galaxies significant 
non-axisymmetric forces were detected in the outer parts of the disks connected to spiral arms, corresponding even to bar 
strength class 3. 

(4) We confirm the previous result by BB and Block et al. (2001) showing the large overlap in bar strength between the 
optical SA, SAB and SB classes. Actually, SB-galaxy can belong to any of the bar strength classes between 1 and 6. 

(5) We found that 95 % of SB galaxies in our sample belong to the "classical bars" identified in the neax-IR, which means 
that the bars are similar. In the optical 1/3 of the "classical bars" are not classified as SB. Even bars which are in the optical 
obscured by dust and which become dust penetrated in the near-IR, cover all bar classes from 1 to 6, thus indicating that 
even strong bars can be obscured by dust. 

(6) Bar lengths are estimated from the phases of the m=2 and m=4 Fourier components of density requiring that the 
phase is maintained nearly constant in the bar region. Bar length is found to correlate with the galaxy brightness Mb, 
confirming the previous result by Kormendy (1979) in the optical region. 

(7) The number of SB-galaxies in our sample drops rapidly at inclinations > 50°. A similar bias appears also in RC3 
when limiting to galaxies brighter than 15.5 uiag. This bias might have important implications while studying the frequencies 
of bars at low and high redshifts, especially because the bias increases toward fainter galaxies. Also, at high inclinations the 
number of SAB galaxies in RC3 drops much more rapidly than the number of SB-galaxies, thus being a manifestation of an 
ambiguous definition of the de Vaucouleurs's SAB class. 
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FIGURE CAPTIONS 

Figure 1. Different models for the vertical density profile are compared. In the upper row the commonly used exponential, 

sech^ and sech laws are compared, with the scale-factor h chosen in a manner that yields identical slopes for z >> hz- The 
frame in the left displays the vertical profiles, while in the right the resulting convolution functions g are shown. These models 
correspond to those applied in Table 4 for calculation of Qb- The lower row compares 5 different model profiles, including 
also Gaussian and uniform models. The scale-factors arc chosen in a manner that yields the same vertical dispersion in each 
case. For comparison, the convolution function corresponding to 2D softened gravity is also shown, with softening parameter 
e = V< 2^ >. These same models are applied for calculation of QT-profiles for NGC 1433 in Fig. 3. 

Figure 2. Testing of the calculation method. The left hand frames show the results of our polar method applied to 
high quality H-band image of NGC 1433 (Buta et al. 2001), using different maximal numbers of even Fourier components as 
indicated in the frame. Fourier components were calculated with two pixel wide zones (1 pix =1.141 arcsec), and the image 
was divided to 128 azimuthal bins. In the upper frame hr/hz = 2.5, typical for early-type galaxies, while the lower frame 
shows the influence of assuming 4 times thinner disk (For hr = 45" and distance of 11.6 Mpc the studied hz values are about 
Ikpc and 0.25 kpc, respectively). For comparison, the right-hand frames show results using a Cartesian evaluation of the 
gravitational potential, with different number of grid divisions. The results for 512 x 512 Cartesian grids are also super-posed 
on the left-hand curves. 

Figure 3. The effect of using different vertical density laws on the obtained Qr-profiles for NGC 1433. The density 

laws correspond to those in the lower row of Fig. 1 (scale-factors chosen to yield the same vertical dispersion), and the plot 
indicates how progressively more peaked vertical density profiles yield stronger non-axisymmetric forces. The two studied 
values of hr/hz are the same as in Fig. 2.: for thinner disk the precise form of the density profile is less significant. 

Figure 4. The effect of radius-dependent scale-height on Qrir) profiles for NGC 1433. As in Fig 2, an exponential vertical 
profile is assumed for each distance, but with either positive {dhz/dr = 0.05, solid curves) or negative {dhz/dr = —0.05, dashed 
curves) gradient of the vertical scale height with distance. Three different cases are studied, where hr/hz = 2.5, either at 
r = 1,2, or 3hr. The lines indicate the vertical scale height-profile in each case, with hz multiplied by 0.01. According to 
Schwarzkopf & Dettmar (2001), the disk thickness often increases with radius, and dhz/dr = 0.05 corresponds to the maximum 
value observed for early types. Negative gradients are also studied, to emphasize the smallness of the expected maximal effect 
of vertically dependent hz - For comparison, the crosses indicate the peak of the Qt profile in the case of constant hr/hz = 2.5. 
The polar method is used. 

Figure 5. A more detailed comparison of the effects of distance-dependent hz, corresponding to the middle frame in Fig. 

4, except that twice larger positive and negative gradients arc studied, together with the case of constant hz. The mean radial 
force profile is shown, together with the m = 2 and m = 4 Fourier amplitudes of the tangential force components (forces 
are in arbitrary units). As expected, tangential force components are increased in regions where scale height is reduced. The 
same is true also for the radial force, which however is even more strongly affected in the region of maximal Ft/Fr ratio, 
explaining the slightly reduced Qt in the case of positive dhz/dr gradient. 

Figure 6. Bar strength measurements by BB and by us are compared for 13 galaxies in H-band. In both measurements 
the same orientation parameters were used and the vertical mass distribution was approximated by an exponential profile 
with hz = 325 kpc {Ho = 75 km Mpc~^). The m = 0, 2, 4, 6 Fourier components of density were included, and calculated 
for radial annulae with width of 2". In the azimuthal direction 128 divisions were used in all cases. No softening was used in 
our potential evaluation. 

Figure 7. For the "classical barred" galaxies in the near-IR we show: the original 2MASS image in the plane of the sky 
(1*'' column), the same image de-projected to face-on orientation, with the m=0 Fourier component being subtracted (2*'' 
column), the "butterfly pattern" (3*'' column), the m=0,2,4 surface density proflles (4"' column) and the radial QT-profllc 
(5*'' column). The "butterfly patterns" were calculated from the potentials using all even and odd components up to m=6, 
while in the Qr-profiles all even components up to 10 were included. In the butterfly diagrams the thick contours correspond 
to Qt levels of 0.1, 0.2, 0.3, 0.4, 0.5, while the thin contours refer to corresponding negative levels. In the Qr-profiles the thin 
lines show measurements in the four image quadrants, whereas the thick lines show the mean values. The dashed vertical line 
in the Qr-profile shows the length of the bar estimated from the phases of m=2 and m=4 density components as explained in 
the text. The numbers in the upper right corner give the values for the maximum Qt in the bar region and its radial distance. 
The maximum is also shown by the box symbol in the Qr-profiles. In the de-projected image the solid circle indicates the 
location of the maximum Qt and the dashed circle the measurement region. 

Figure 8. Absolute blue magnitude v.s. bar length. The magnitudes are from RC3 and bar lengths are estimated from 
the phases of the m=2 Fourier amplitudes as explained in the text. In Fig. 8a bar lengths are given in absolute units, using 
the distances from TuUy (1988), and in Fig. 8b they are scaled to the scale-length of the disk. H = 75 km s~^Mpc~^. 

Figure 9. The inclination distribution of the galaxies in our sample, SB galaxies and the "classical bars" identified in 
the near-IR (see the text) being shown separately. For comparison, the inclination distribution for all galaxies in our sample 
is also shown. The dashed line indicates a random distribution of inclinations, normalized to the maximum inclination and 
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the number of galaxies in the sample. 

Figure 10. A correlation between the inclination of the disk (INC) and bar strength (Qb) for the galaxies in our sample. 

Figure 11. A similar figure as Fig. 9, but shown for the de Vaucouleurs's (1963) classes SA, SAB and SB for a magnitude- 
limited sample (Bt < 15.5 mag) from RC3. 

Figure 12. Distribution of the non-axisymmetric forces among the de Vaucouleurs's (1963) classes SA, SAB and SB, 
measured in J, H and K-bands (Figs a, b and c, respectively). For comparison, bar strength distribution is shown also for the 
"classical bars" identified in the near-IR. Notice that not all non-axisymmetric forces in the figure are associated to bar-like 
potentials, being due to spiral arms for some SA and SAB-galaxies. 
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Table 3. Bar strengths for galaxies without classical bars. 
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S = Qt due to spiral arms B = bar-like "butterfly" structure, but no clear Qr-maximum in the assumed bar 
region BM= both Q-p-maximum and bar-like butterfly structure appear 
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Table 4. For the test 


galaxies: bar strengths and the parameters 


used in the calculation. 
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